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Abstract 

The thermodynamic properties of (TeFg) 59 clusters that undergo temperature-driven phase transitions have been 
calculated with a canonical J-walking Monte Carlo technique. A parallel code for simulations has been developed 
and optimized on SUN3500 and CRAY-T3E computers. The Lindemann criterion shows that the clusters transform 
from liquid to solid and then from one solid structure to another in the temperature region 60-130 K. 
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1. Introduction 



tential: 



The interest in phase transitions in finite systems 
is two-fold: firstly, bulk properties of the material 
can be simulated if the system is studied under pe- 
riodic boundary conditions [1] and proper account 
for rounding and shifting of the measurable quan- 
tities is taken in the analysis. Secondly, free finite 
systems, such as molecular clusters, are of interest 
due to their peculiar properties, which are not ob- 
served in bulk systems of the same substance [2]. 
These are important for the synthesis of nanoma- 
terials and nanodevices [3] . 

We use a canonical Monte Carlo (MC) method to 
investigate temperature-induced critical behavior 
in clusters containing 59 TeF§ molecules. 

Our model system contains 59 octahedral rigid 
TeF 6 molecules which are allowed to rotate and 
translate. Their interaction is reliably described by 
a Lennard-Jones and a Coulomb atom-atom po- 
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where q is a generalized coordinate; is the dis- 
tance between the i — th and j — th atom. The 
indices a, (3 denote cither a fluorine or a tellurium 
atom. The parameter values are taken from [4] . 

This orientation-dependent potential generates 
rugged potential energy surfaces. How to tackle 
the behavior of a system having a rugged potential 
energy surface in MD simulations is discussed in 
[5,6]. 

Here we tackle the problem with the help of the 
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jump- walking Monte Carlo method, [7,8], which is 
a smart development of the Metropolis algorithm. 
We have parallelized the code to speed up the com- 
putations [9]. 

2. J- walking algorithm 

In its general form, the method generates trial 
moves from a higher-temperature (Xj) equilib- 
rium distribution with a probability specified by 
the variable Pj: 

P = min{l,q[x*\x]} (2) 
where the Metropolis probability function is: 

q(x*\x)=e-(±-^ V ^- V ^ (3) 

V(x) is the potential energy of the system at a 
temperature T and V(x*) is the potential energy 
of a new configuration at the temperature Tj. 

The remaining trial moves are of conventional 
Metropolis character. 

In our study we use multi-temperature general- 
izations [10] of the basic approach. 

We develop a parallel J-walking code that en- 
ables us to carry out a Monte Carlo simulation 
efficiently in a multiprocessor computing environ- 
ment. We apply the code in the study of the ther- 
modynamic properties of (TeF^n, n = 59 clusters. 

The J-walking technique can be implemented in 
two ways. The first approach is to write the con- 
figurations from the simulation at the J-walking 
temperature to an external file and access these 
configurations randomly, while carrying out a sim- 
ulation at the lower temperature. It is necessary to 
access the external files randomly to avoid corre- 
lation errors [11]. The large storage requirements 
limit the application of the method to small sys- 
tems. The second approach uses tandem walkers, 
one at a high temperature where Metropolis sam- 
pling is ergodic and multiple walkers at lower tem- 
peratures. 

The best features of these two approaches can 
be combined into a single J-walking algorithm [12] 
with the use of multiple processors and the Mes- 
sage Passing Interface (MPI) library. We incorpo- 
rate MPI functions into the MC code to send and 



receive configuration geometries and potential en- 
ergies of the clusters. Instead of generating exter- 
nal distributions and storing them before the ac- 
tual simulation, we generate the required distri- 
butions during simulation and pass them to the 
lower-temperature walkers. 
Parallel J-walking algorithm: 

- Step 1. For each t make Metropolis MC steps: 

• Rotate each molecule. 

• Translate each molecule. 

• Reject or accept step. 

• Go to step 1. 

- Step 2. After Si steps - collect statistics: 

• Potential energy histogram. 

• Energy average and deviation. 

• Heat capacity Cy. 

• Save current configuration. 

- Step 3. After S2 steps - make jump- walking step 
by exchanging the configurations using MPI. 

- Step 4. Go to step 1. 

The diagram of model J-walking is shown in fig. 1 . 
Each square is a Metropolis MC simulation at 
a particular temperature. The set of boxes on 
the right-hand side represents the array of pre- 
vious configurations of the system, which are 
stored in the memory to avoid correlations be- 
tween the lower and the higher temperature. At 
each trial jump we randomly choose one of the 
4 systems. When a configuration is transmitted 
to a lower-temperature process, it is a configura- 
tion randomly chosen from the array of higher- 
temperature walkers. The current configuration 
of the walker then replaces the configuration just 
passed from the array to another temperature. 
In fig. 2 we show the parallel decomposition of 
computation. Each process computes part of one 
of 4 multistage J-walking chains, and exchanges 
configurations and energy with others. 

We use array sizes of 2500 configurations. The 
number of configurations is limited by the proces- 
sor RAM used in simulation. In the parallelization 
implemented in our code, the arrays are small and 
do not inhibit applications of the method to large 
systems. For the computations in the present work 
the number of MC passes for each temperature is 
6.10 5 for each cluster containing 59 TeF§ molecules 
(413 atoms). We make J-walk jump attempts at 
every 50 Metropolis MC steps during the thermali- 
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sation and at every 150 steps during the main com- 
putation. 

The computer code has been ported, tested and 
optimized on SUN 3500 and CRAY-T3E machines. 
In our program dynamic memory management has 
been implemented for optimal usage of the mem- 
ory. We find that memory and performance re- 
quirement make CRAY-T3E more suitable for such 
computations. In our runs we use 64 processors 
each with 64 MB RAM. Each run of 6.10 5 steps 
takes approximately llh per CPU. 



3. Results and Conclusions 

Using the parallel code described in the previ- 
ous section we make sets of production runs for 59 
molecule clusters. The previous MD analysis of the 
temperature behavior of TeF§ clusters pointed out 
a two-step structural transformation process [13] 
from an oricntationally disordered bcc structure 
below cluster solidification to an oricntationally 
oriented ordered monoclinic structure detected at 
low T (below 20K). The first step involves lattice 
reconstruction (bcc to monoclinic) and a partial 
order of one of the molecular axes, when the cluster 
is cooled down to its freezing point. This transition 
has been proved to be a first-order phase change 
[14] by detecting coexisting phases. A further tem- 
perature decrease causes complete orientational or- 
der of the three molecular axes. This transforma- 
tion is continuous. The diagnostic method devel- 
oped in [13] has been implemented to animate the 
solid-solid transformations [15]. 

To distinguish between the different phases, 
which the clusters adopt at different temperatures, 
we have computed the Lindemann index [16]: 
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where rij is the distance between the centers of the 
molecules (Te atoms). For the values of Su n > 0.08 
the system is in liquid phase and for 5u n < 0.08 is in 
the solid phase. Fig. 3 presents the Lindemann cri- 
terion for the 4 systems. The solidification occurs 



in interval 100if -llbK . In the interval 75K-80K 
the angle of the curve has changed. 

We calculate the radial distribution function to 
reconstruct the crystal lattice. Fig. 4 and fig. 5 show 
the normalized radial distribution function for dis- 
tances between centers of the TeFg molecules (Te 
atoms). Fig. 4 is at 157. The peaks of the curve 
correspond to disordered liquid structure. Fig. 5 is 
at 76K where the peaks correspond to an ordered 
structure (bcc lattice). 

In this study a J-walking Monte Carlo algorithm 
has been implemented to study (TeF^) clusters. A 
parallel J-walking code has been developed. Ther- 
modynamic properties have been observed in the 
range from 60.fr to 150-fT. We can compare the re- 
sults to the results obtained from MD simulations. 
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Fig. 1. Diagram of a simple J- walking process. T\ and T2 
present the systems at two different temperatures 



□ [mid i 





* L * 



□ [Tumi n nnmn 

i l _ 



□ nTTTTTTl 




T, : □ imTTTTI ! i □ MM 



T 4 1 □ nTTTTTTl ; ; □ TJHUJE : ; □ MM 



□ nTrjrjTj \ 



Fig. 2. Diagram of a parallel decomposition of data and 
calculation between processors in J-walking algorithm. In 
this example we have 3 different systems at temperatures 
Ti , T2, T3 and T4 distributed between 6 processors 
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Fig. 3. Lindemann criterion calculated for 4 clusters of 59 
molecules. 
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Fig. 4. Radial distribution function for 59 molecule TeFg 
cluster at 157AT. 
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Fig. 5. Radial distribution function for 59 molecule TeF@ 
cluster at 76 K. 
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